Critical behavior of frustrated spin systems with nonplanar orderings 
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The critical behavior of frustrated spin systems with nonplanar orderings is analyzed by a six-loop 
study in fixed dimension of an effective 0(N)xO(M) Landau- Ginzburg- Wilson Hamiltonian. For 
this purpose the large-order behavior of the field theoretical expansion is determined. No stable fixed 
point is found in the physically interesting case of M = N = 3, suggesting a first-order transition 
in this system. The large N behavior is analyzed for M = 3, 4, 5 and the line N c (d = 3, M) which 
limits the region of second-order phase transition is computed. 
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I. INTRODUCTION 

The critical behavior of frustrated spin systems with 
noncollinear order represents an important and, at the 
same time, enchanting issue d' 213 The debated topic con- 
cerns the nature of the phase transition which these 
systems should undergo. In the case of a second-order 
phase transition, the universality class is expected to be 
noveli*^ and characterized by new critical exponents dif- 
ferent from those of 0(7V)-symmetric models. 

Typical systems which have been studied for many 
years and where the situation is still unclear are the 
staked triangular antiferromagnets and helimagnets, 



which present a noncollinear but planar spin orderings 
of the ground statei^ On this issue, there is much de- 
bate since field theoretical methods, Monte Carlo sim- 
ulations and experiments provide contradictory results 
in many cases, sometimes supporting a weak first-order 
phase transition, and sometimes a second-order one (see 
Refs. IH3 for reviews on this topic). The perturbative 
field theoretical studies on canted magnetic systems are 
based on an effective 0(N)xO(M) Landau-Ginzburg- 
Wilson (LGW) Hamiltonian with M < N, which de- 
scribes M-dimensional spin orderings in isotropic iV-spin 
spac o 1 : 2 : 4 with the O (N)/0 (N — M) pattern of symmetry 
breaking: 



J 



where <f> a (1 < a < M) are M sets of iV-component vec- 
tors. Negative values of vq correspond to simple ferro- 
magnetic or antiferromagnetic ordering, and to magnets 
with sinusoidal spin structures i The condition < vq < 
M I (M — l)ito is required to have noncollinearity and the 
boundedness of the free energy. For M = 2 the Hamilto- 
nian m describes magnets with noncollinear but planar 
orderings as the XY (N — 2) and Heisenberg (N = 3) 
frustrated antiferromagnets, while for M > 3 systems 
with noncoplanar ground states i** 

The physical relevant case M = N = 3 represents the 
noncoplanar spin orderings which are three-dimensional 
in spin space and can be realized in structures in which 
the instability occurs simultaneously at three inequiva- 
lent points in the wavevector space as shown in Fig. 

For M = 2, the critical properties of the Hamiltonian 
JTJ have been extensively studied in the framework of the 
e = 4 - d^&L l/jV ) 1 ^?!? I = d - 2 expansions^ as well 
as in fixed dimension d = 



, and d = 3 t 1 'V 4 i 15 i 16 i 1 7 
and with the exact renormalization group (ERG) tech- 
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niqueiiSiiSiSSi The existence and stability of the fixed 
points were found to depend on N in three dimensions. 
For sufficiently large values of components of the order 
parameter four fixed points exist in the region u, v > 0: 
Gaussian, 0(2 AT), chiral and antichiral. The chiral fixed 
point is the only stable one. Decreasing the value of N 
the e expansion predicts the existence of an N C (M = 2) 
greater than three at which the chiral and antichiral fixed 
points coalesce, disappearing for TV < N C (M = 2). These 
results, which are in agreement with the and 
the ERG oneSfi^ suggest the existence of a first-order 
phase transition for XY and Heisenberg frustrated an- 
tiferromagnets. However one should consider that both 
the e and 1/N results are essentially adiabatic moving 
from large ./V and small e. For this reason they do not 
necessarily provide the (essentially nonperturbative) fea- 
tures of the models in the region below N C (M — 2). In 
fact six-loop calculations performed directly in three di- 
mensionsi 4 *^ showed a different scenario characterized 
by two critical numbers of components of the order pa- 
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FIG. 1: Examples of nonplanar spin orderings (M — 3) in fee 
antiferromagnets of type-II (left) and type-I (right). 



rameter N C (M = 2) = 6.4(4) and N c2 (M = 2) = 5.7(3): 
for N c2 {M = 2) < N < N C (M = 2) no stable fixed point 
is present in the RG flow, while for TV < N C 2(M = 2) a 
stable chiral fixed point appears again^ which is a fo- 
cus^ not related with its counterpart found within the e 
and 1/N expansions. The discrepancy between the ERG 
predictions and the weak-coupling ones remains an open 
problem which could be explained either with the low 
order considered in the derivative expansion in the first 
case, or by the problem of Borel-summability in the sec- 
ond one. 

For generic values of M > 3, the e and 1/N expansions 
give rise to the same scenario with the only difference 
being the value of N C (M) which varies with M. In the 
physically important case of M = N = 3, the e expansion 
predicts the absence of a stable chiral fixed point. Sev- 
eral Monte Carlo simulations2ii22i2ii2i based essentially 
on Stiefel's models and ERG results^ agree indicating a 
first-order phase transition for this system. Anyway no 
data in the fixed dimension expansion exist. 

In this paper we analyze the Hamiltonian |JTJ for M > 3 
through a six-loop perturbative study in the framework 
of the fixed dimension approach, in order to investigate 
the critical behavior of frustrated spin systems with non- 
planar orderings directly in three dimensions and to have 
a global picture of the region where the results obtained 
with different perturbative methods are in agreement. 

The paper is organized as follows. In Sec. [HI we in- 
troduce the field-theoretical approach for the effective 
LGW Hamiltonian (JTJ, in Sec. IIIII we present the analy- 
sis method and we discuss the singularities of the Borel 
transform, in Sec. II VI we present the results and finally 
in Sec. 0we draw our conclusions. 



II. FIXED DIMENSION PERTURBATIVE 
EXPANSION OF THE 0(N)xO(M) LGW 
HAMILTONIAN 

The Hamiltonian is studied in the framework of 
the fixed dimension field-theoretical approach. As in the 
case of XY and Hcisenberg frustrated models^, the idea 
is to perform an expansion in terms of the two quartic 
coupling constants uq and vq, and to renormalize the 
theory imposing the following conditions on the (one- 
particle irreducible) two-point, four-point and two-point 



with an insertion of ^(f> 2 correlation functions: 



lfi» = S^SijZ^ 1 [m 2 +p 2 + 0(p 4 )} , (2) 



^i,i,bj,ck,dl (P) — ~zSai,bj,ck,<U + T C a i,bj,ck,dl , (3) 



(4) 



where S al ,b 3 ,ck,di an d C al ^,ck,di are appropriate combi- 
natorial factorsiii*^ 

The perturbative knowledge of the functions r^, I^ 4 ) 
and r^ 1 ' 2 ) allows one to relate the renormalized parame- 
ters (u,v,m) to the bare ones (uq,Vq,t ). 

The fixed points of the model are defined by the com- 
mon zeros of the (3 functions, 



/3 u (u, v) =m 



du 
dm 



(3 v (u, v) = m 



dv 
dm 



(5) 



and the stability properties of these points are deter- 
mined by the eigenvalues u>i of the matrix: 



/ d(3 u (u,v) df3 u (u,v) \ 
du dv 



\ 



d/3 v {u,v) d/3 v (u,v) 



du 



dv 



(6) 



/ 



A fixed point is stable if both the eigenvalues have a pos- 
itive real part, while a fixed point which possess eigenval- 
ues with nonvanishing imaginary parts is called a focus. 

The values of the critical exponents 77, v are related to 
the RG functions 77^, and rj t evaluated at the stable fixed 
point: 



V = V<t>{u ,v ), 
v = [2 - 77^ (it*, v*) + rj t (u*, v*)] 1 



where 
Vt(u,v) 
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du 

dlnZt 
du 



■0v 



(7) 
(8) 

d_^ m 

dlnZ t 

dv 



(10) 



The other critical exponents can be derived by scaling 
relations from 77 and v. 

In this work the following symmetric rescaling of the 
coupling constants is adopted in order to obtain finite 
fixed point values in the limit of infinite components of 
the order parameter: 



_ 167T _ 

u = -—-Rmnu, 



167T 

v = -—Rmnv (11) 



where R N = 9/(8 + AT). 

For N — ► oo four fixed points exist: Gauss (u = 0, v = 
0), Heisenberg (w = l,u = 0), chiral (u = M, v = M) 
which is stable, and antichiral (u = M — l,v = M). In 
the following we denote with and /3 C the (3 functions 
corresponding to the above rescaled couplings. 



III. RESUMMATION AND ANALYSIS 
METHOD 

The perturbative six-loop series^ for the RG func- 
tions are asymptotic and some resummation procedure 
is necessary in order to extract quantitative physical in- 
formations. In this work the property of Borel summa- 
bility is exploited resumming the perturbative series by 
a Borel transformation combined with a conformal map- 
ping^ (CM) which maps the domain of analyticity of the 
Borel transfor m (cut at the instanton singularity) onto a 
circle (see Refs. l26l27l for details). This resummation pro- 
cedure needs the knowledge of the large order behavior of 
the perturbative series which is connected with the sin- 
gularity (flb) closest to the origin at fixed z = v/iii 11 ^ 5 
For the Hamiltonian Q with generic 0(N)xO(M) sym- 
metry and M < N we find 



1 r 2M 

— = -ciRmn for < z < — — -, 
Uf, M — 1 

— = —ciRmn 

Ub 
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forz < 0, z > 



(12) 

2M 
M — 1 ' 



where a = 0.147777422. 

In the case of M = 2 the singularity reduces correctly 
to the one of Ref. 0. For z > j^zj there is a singularity 
on the real positive axis which however is not the closest 
one to the origin for z < jP^r . Thus, for z > jMrr the 
series are not Borel summable. The region where the 
perturbative series are not Borel summable is the same 
of the one where the condition for the boundedness of 
the free energy is satisfied. It is easy to show that the 
conditions are equivalent. 

With this resummation procedure one obtains a lot of 
approximants for each RG perturbative series character- 
ized by three parameters p, b and a which can be varied 
in order to estimate the systematic errors in the final 
results. If R(u,v) is a perturbative series in u and v 



l l-k 



R(u,v) = J^iW 1 ^, 

fc=0 h=Q 

the approximants will have the following form: 

p 

E(R) p (a, b; u,v) = B k (a, 6; v/u) 



(13) 



fc=0 



y(ut; v/u) 




FIG. 2: Zeroes of the /3-functions for N = 3 and M = 3 in the 
(u, v) plane. Pluses (+) and crosses (x) correspond to zeroes 
of /3b (u, v) and f3u(u,v) respectively. 



where 



y{u;z) 



y/1 - u/ub{z) - 1 

y/l - U/Tl b (z) + 1 



(15) 



Jo [l-y{ut-v/u)Y 



(14) 



and the coefficients Bk arc determined by the condition 
that the expansion of E(R) v (a, b;u,v) in powers of u and 
v gives R(u, v) to order p. 

In order to find the fixed points of the Hamiltonian 
pjl the two (3 functions are resummcd in the whole zone 
(u > 0, v > 0). For each (3 function the approximants 
are chosen which stabilize the series for small values of 
the coupling constants with varying the considered per- 
turbative order (number of loops), i. e. a = 0,2,4 and 
b = 3,6,9,12,15,18. The stability eigenvalues are 
evaluated by taking the approximants for each (3 func- 
tion and by computing the numerical derivatives of each 
couple of them at their common zero. Only the approx- 
imants which yield fixed point coordinates compatible 
with their final values are considered. The critical expo- 
nents are computed by choosing the approximants which 
are more s table with varying the perturbative order (see 
Ref. Illll2l for details). 



IV. RESULTS 

Fixing M = 3, the analysis of the six- loop RG se- 
ries reveals the absence of a stable chiral fixed point for 
the physical interesting case of N = 3; in this case no 
second-order phase transition is expected as it is clear 
from the picture of the zeros of the (3 functions in Fig. 
121 In this Figure the domain 0<u<6,0<u<6is 
divided in 40 2 rectangles, and all the sites in which at 
least two approximants for (3 a and (3y vanish are marked. 
These results are in agreement with the ERG prediction 
of Ref. 19 where no stable fixed point is found, and with 
Monte Carlo simulations on the corresponding Stiefel's 
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mode! 23,2 which show a first-order phase transition. In- 
creasing the value of components of the order parameter 
N, the curves of the zeros of the two (3 functions become 
closer and it exists a value N C (M = 3) at which they in- 
tersect each other. For N > N C (M = 3) four fixed points 
exist in the region v > 0,u > 0; the Gaussian one, the 
Heisenberg one, the antichiral one and the stable chiral 
one. This last drives the system to a continuous phase 
transition. For N = N C (M = 3) the chiral and antichi- 
ral fixed point merge and for N < N C (M = 3) no stable 
fixed point is found. The estimate of iV c is obtained by 
considering all the 324 possible combinations of the ap- 
proximants for the two (3 functions and evaluating for 
each of them the value of N c . With this procedure we 
have N C (M = 3) = 11.1(6) at six-loop order. 

The same scenario of fixed points holds for M = 4 and 
M = 5 with different values of N c . 

All the five-loop and six-loop estimates of N c for 
M = 3, 4 and 5 are shown in Tab. [I] where the e ex- 
pansion and the 1/N predictions of Ref. ?j are reported 
for a comparison. The agreement of the results is good 
even if it worsens with increasing the magnitude of M. 
Nevertheless, in the case of a comparison with the e ex- 
pansion results, only three-loop series are known which 
need to be resummed in order to extract some informa- 
tion. Due to the low perturbative order considered, one 
may only have an indication of the right results. In fact 
the perturbative series can be evaluated extending them 
analytically by means of a Pade Borel transform^ (e. g. 
N C (M = 2) ~ 3.39) or by different techniques^ (e. g. 
N C (M = 2) ~ 5.3) and the results strongly depens on the 
method applied. In the case of the 1/N expansion, one 
has to consider that N c is determined by perturbative se- 
ries which are known only to 0(1/N 2 ), the correction to 
the leading term is not small, and since N c is expected not 
to be large in magnitude, one is extrapolating the results 
to value of N which could be dangerous. Furthermore 
the results of the different perturbative approaches are 
nearer in magnitude for small values of M because with 
increasing the perturbative order considered, the correc- 
tions to the leading terms become more important for 
bigger M. Even the fixed-dimension results become less 
stable with increasing the value of M. In order to provide 
a complete characterization of the critical properties of 
the Hamiltonian Q we analyze it for large values of com- 
ponents of the order parameter at fixed M. The location 
of the chiral fixed point, the stability eigenvalues and the 



TABLE I: Comparison between values of N a (d = 3, M) ob- 
tained with different perturbative methods. 



M 


2 


3 


4 


5 


CM 5-loop 




11.6(6) 


15.4(9) 


19.0(1.2) 


CM 6-loop 


6.4(4)i4 


11.1(6) 


14.7(8) 


18(1) 


e-expansion— 


5.3(2) 


9.1(9) 


12(1) 


14.7(1.6) 


1 /iV-expansion 7 


5.3 


7.3 


9.2 


11.1 



critical exponents are calculated for TV = 32, 64 and are 
displayed in Tab.lTTl For the evaluation of the critical ex- 
ponents 7 and v we resum the perturbative series of 1 /'y 
and 1/v because of their better behavior. In these cases 
one can see a quite good agreement or the results with 
the 1/N expansion predictions. The data of Tab. [Q] cor- 
roborate the expectation for which the line N C (M) limits 
the region where all the perturbative investigations bring 
to the same predictions. 



V. CONCLUSIONS 



In this work we have studied the critical thermody- 
namics of the 0(N)xO(M) spin model for N > M > 3 
and M < 5 on the basis of six-loop RG series in fixed 
dimension. We have evaluated the large order behav- 
ior of the perturbative RG series and we have analyzed 
them by means of a CM resummation technique. For 
the physically interesting case of N = M = 3, no stable 
fixed point exists in the RG flow diagram leading to the 
conclusion that this system is expected to undergo a first- 
order phase transition. This result agrees with previous 
ERG predictions^ and Monte Carlo studies on Stiefel's 
modelsiS 3 - 2 ^ For M = 3, 4, 5 we have computed the crit- 
ical value N C (M) of components of the order parameter 
which separates the region of first-order phase transition 
to the second-order one, finding N C (M = 3) = 11.1(6), 
N C (M = 4) = 14.7(8) and N C (M = 5) = 18(1). In 
this case the comparison with the and 1/N expan- 
sioni predictions is quite good even if it worsens with 
increasing the magnitude of M. Finally, being the region 
N > N C (M) expected to be connected with that accessi- 
ble by the large N analysis, we have given a characteri- 
zation of the critical properties of these systems for large 
values of N finding a satisfactory agreement between the 
two perturbative methods. 
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TABLE II: Large- N results. 



AT = 64, M = 3 


Ch. 


[u>i, W2] 


V 


7 


CM. (6-loop) 


[2.934(1),3.0051(8)] 


[0.960(1),0.915(2)] 


0.960(2) 


1.911(4) 


CM. (5-loop) 


[2.935(2) ,3.006(2)] 


[0.958(2),0.911(3)] 






i / i\ t • 7 8 

1 /A'-expansiori«- 




[0.97,0.93] 


0.965 


1.92 


N = 32, M = 3 


Ch. 




V 


7 


CM. (6-loop) 


[2.861(2),3.012(2)] 


[0.919(2),0.810(6)] 


0.912(4) 


1.807(8) 


CM. (5-loop) 


[2.862(7),3.014(7)] 


[0.916(5),0.80(1)] 






1 / AT • 7 8 

1 /A'-expansioi: 




[0.94,0.86] 


0.927 


1.84 


A = 64, M = 4 


Ch. 




V 


7 


CM. (6-loop) 


[3.8550(15),3.950(1)] 


[0.9536(11),0.891(2)] 


0.949(3) 


1.886(6) 


CM. (5-loop) 


[3.856(4),3.951(3)] 


[0.952(3),0.886(4)] 






i / i\ r • 7 8 

1/iV -expansion ' 




[0.97,0.91] 


0.956 


1.90 


N = 32, M = 4 


Ch. 


[wi, W2] 


V 


7 


CM. (6-loop) 


[3.7035(30),3.905(3)] 


[0.905(2),0.75(1)] 


0.887(6) 


1.75(1) 


CM. (5-loop) 


[3.65(7) ,3.906(11)] 


[0.904(6)0.73(2)] 






i / ?vr • 7 8 

1/iV -expansion 




[0.93,0.83] 


0.907 


1.80 


AT = 64, A/ = 5 


Ch. 




V 


7 


CM. (6-loop) 


[4. 7676(16) ,4.8864(15)] 


[0.9467(12), 0.800(3)] 


0.937(3) 


1.862(6) 


CM. (5-loop) 


[4.770(6) ,4.888(5)] 


[0.945(3),0.860(7)] 






1 /W-expansior. 7,8 




[0.96,0.90] 


0.946 


1.88 


N = 32, M = 5 


Ch. 


[Wl, W2] 


V 


7 


CM. (6-loop) 


[4.528(6) ,4.783(3)] 


[0.8917(20),0.675(15)] 


0.860(8) 


1.696(16) 


CM. (5-loop) 


[4.52(2),4.78(2)] 


[0.893(6),0.65(2) 






1 /W-expansion-i- 




[0.92,0.80] 


0.887 


1.75 
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